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We study the target normal spin asymmetry for the reaction ep — > ep at electron energy up to few 
GeV in the laboratory frame. The asymmetry is proportional to the imaginary part of the reaction 
scattering amplitude. To estimate the imaginary part of the amplitude we use the unitarity relation 
and saturate the intermediate hadron states by the proton and resonances from the first, second and 
third resonance regions. The resonance electromagnetic transition amplitudes, needed to evaluate 
the asymmetry, are taken from experiment. 
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I. INTRODUCTION 

A study of the elastic ep-scattering is an important source of information about the internal structure of the proton. 
Due to the smallness of the fine structure constant, a w 1/137, the first order perturbation term (the one-photon 
exchange) is assumed to give a main contribution to the electromagnetic transition amplitude. In the one-photon 
approximation the elastic ep-scattering is described by two quantities, the electric, Ge(Q 2 )-, and magnetic, Gm(Q 2 ), 
form factors. 



> 

' The form factors Ge and Gm are usually extracted from cross section data by the Roscnbluth separation method, 
f*"*- . The database for Ge and Gm obtained by this method shows that the ratio Ge/Gm is approximately constant. 

Recently new precise measurements of the ratio Ge/Gm were done at Jefferson Lab PHU Q by the recoil polarization 
*n , method 0, 0] and yield significantly different results than the Rosenbluth separation method . 

Because both measurements are based on the one-photon exchange approximation it is natural to assume that this 
discrepancy may be explained by the second order perturbation term in the a expansion. The a 2 perturbation term 
should result in many effects. Its real part contributes to the cross-section and destroys the Roscnbluth formula. The 
imaginary part appears in one-particle polarization observables of the elastic ep-scattering, the target and beam spin 
asymmetry, which vanish in the one-photon approximation. The present experimental technique makes possible to 
' measure such observables and thus to take under control effects beyond the one-photon approximation. 

The aim of the present paper is a calculation of the target spin normal asymmetry in the elastic ep-scattering at 
Ch ■ electron lab. energy E\ a h < 2 GeV. 

The imaginary part of the scattering amplitude, which determines the asymmetry, is simply related (through 
the unitarity condition) to the electroproduction amplitudes of different hadronic states. The so-called "elastic" 
contribution (i.e. when the hadronic intermediate state, entering the unitarity condition, is the proton) to the 
asymmetry was calculated in ■ In Q authors got strict bounds on the "inelastic" part of the asymmetry using the 
Schwartz inequality. But, as the authors noted themselves, those bounds highly overestimate the actual values of the 
asymmetry, especially at high scattering angles. 

In the recent work [(| the asymmetry was calculated with N and nN intermediate states. Such approach gives a 
reasonable approximation at low electron energies, but becomes worse as the energy increases. 

Contrary to ||, we will calculate the contribution of the resonances in the intermediate state, namely, ^33(1232), 
-Di3(1520), 5*11(1535), Fi5(1680) and Pn(1440), using their experimental electroproduction amplitudes. Such ap- 
proach may be justified at intermediate energy, £i a b ^2 GeV. For the large momentum transfer region parton model 
calculations would be more adequate [Toj . 

It was noted in [j| , that the single-spin asymmetry is sensitive to the electroproduction amplitudes in a wide range 
of photon virtualities, and this may be a new way to access resonance transition form factors. The results of the 
present work may be useful in planning of such experiments. 

The paper in organized as follows. In Section[n]we derive general formalism for the asymmetry, in Sections lIIII and lTvl 
we explain how we describe the resonances and fit their electromagnetic transition amplitudes. Numerical results and 
conclusions are given Section El 
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II. GENERAL FORMULAE 



A. Notation 



We denote the initial electron and proton momenta k and P, respectively, and the final momenta k' and P' . The 
transferred momentum is q = k — k' (q 2 < 0), and the cm. energy squared is s = (P + k) 2 = (P' + k') 2 . Time and 
space components of 4- momenta are denoted like P = (ep, P). M is the proton mass, the electron mass is neglected. 
We denote Dirac matrices 7 M and use the short-hand notation a for . 

Proton spinors with definite helicity A and momentum P are 

, m / V^+M Wx \ , ( e-WcoB g+'(V»-A) \ 



Vep - M (riff) w A /' I e ^/ 2 sin 0+7r( /2 ~ A) 
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where 9, cp are spherical angles of the vector n = P/\P\, a are Pauli matrices. 
Electromagnetic current matrix elements for the proton read 



(P'X'IJ^PX) = UX'(P') U X (P) = Uy{P') 



p 



2M(G e -G m )-± + GmY 

+ 



MP), (2) 



where |PA) is the proton state with momentum P and helicity A, P + = P + P', Ge = GE(q 2 ) and Gm = G]\j(q 2 ) 
are the proton elastic form factors. 

Current matrix elements between the proton (with momentum P) and other hadronic states (with momentum 
P" = P + q) can be expressed via 3 independent invariant amplitudes. In the rest frame of the hadronic state (see 
Figure 0) one has 

£ ( A )(^A|^|PV 2 ) = /f V)<W/ 2 , & X) (h -A|J"|P - Va) = Wf V)«A,A+V a - (3) 

Here % = 7T/ l e l7r< - s ' l_ / 2 ), /i - some hadronic state, s/j, 7r^ - its spin and parity, A - spin projection onto the vector P. 
The quantities can be considered as helicity amplitudes of the process j*p — > h. 

Polarization vectors of a virtual (space-like) photon £ M are defined according to 0. In the coordinate frame of 
Figure n they are 

4°> = _L=(|^| I J > -g°), 4 ±X) = 4(0,Tl,-i,0), (4) 

where the upper index of e shows the spin projection onto the z— axis. 

If the coordinate system is oriented arbitrarily, so that n = —q/\q \ = (cos^sin^, sin<^sin0, cos 6), then 

4 0) = ^y(l<a ~1 Q n), (5) 
eji = {(0,isin(p, — icoscp, 0) =p (0, cos ip cos 6, sin ^3 cos 6, — sin#)}. 
Orthogonality relations 

EW^^v-^f -d rfH A '» = H)V. (6) 
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B. Asymmetry 



The name "target normal asymmetry" corresponds to the situation when the target proton is polarized along the 
normal to the reaction plane, and other particles are unpolarizcd. Under such condition the proton spin has two 
possible directions, say, above and below the reaction plane. Its invariant spin 4-vector should be either collincar or 
anti-collincar with 4-vcctor 



2e v ^ T k v P a P' T 
^q 2 ( S q 2 -(s-M 2 ) 2 ) 



(7) 



which is orthogonal to all momenta and satisfies S 2 = — 1. The corresponding cross-sections, cr-f and crj_, are equal in 
the one-photon approximation, so the difference between them is due to higher order perturbative terms. The target 
normal asymmetry is defined as dimensionlcss ratio 



+ a l 



(8) 



The asymmetry is proportional to the imaginary part of the scattering amplitude. The imaginary part, in turn, can 
be expressed through the unitarity condition, which reads 



(9) 



where i and / are initial and final states, respectively, n is so-called intermediate state and Tfi are T-matrix elements. 
In our case we can, as the first approximation, use one-photon exchange amplitudes in the right-hand side of JHJl. 
Then we obtain 



k k'/ 



2 Im 




S 





0(a s 



(10) 



We denote q% = k — k" q^ = k' — k" and the mass of hadronic intermediate state W = \/P" 2 . 

As it was shown in [8j using (|10fl and the time-reversal symmetry of the electromagnetic interaction, the asymmetry 
can be written as jlfij 



laq 
2n 2 D 



S 

f a _ T/T/ 2 f 1 

- / — — dw 2 / <m k „—L a ^ V w^p'x^px^u^px-^sr^u^p') 

M 2 a p> a p 



in the first non- vanishing order of a, where T a is defined in J3J|, 

' (2s + q 2 - 2M 2 ) 2 



D = 4 



(AM G E - q G M ) + q 2 (AM 2 G% + q 2 G 2 M ) , 



(11) 



(12) 



L""" = Tr(k'^k"^kj a ), 



and the hadronic tensor W^u is defined as 



W^(P'\' p - PX P ) = ^(2^) 4 <S(P + k- P" - k"){P'X' p \J„\h){h\J u \PX p ) 



(13) 



(14) 



Here \h) are all possible hadronic states, which we will refer to as "intermediate states". They can be N, irN, irirN, 
r)N and so on. ^2 is the short-hand notation for 



^ ^ J n (27r) 3 2e, 

N spins' 7 o=l V ; f 



(15) 
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where N is the total number of particles, p a are their momenta. 

Now we will express the hadronic tensor through the electroproduction helicity amplitudes fj^ . 

Consider this tensor in the rest frame of the hadronic state \h), i.e. at P" = 0. We choose coordinate system so 
that both vectors P and P' lie in the yz-plane (see Figure |2J, and the angle between them is f3, < j3 < n. 

The sum E can be split into 2 parts: first, the sum over total angular momentum (=spin) projections, second, 

h 

the sum over all remaining quantum numbers, which we denote E ■ The behaviour of the state \h) with respect to 

h 

spatial rotations is completely described by its spin and spin projection, and does not depend on any other quantum 
numbers. Thanks to this fact, the sum over spin projections can be done explicitly: 

W^P'X'^PXp) = £'E(27r) 4 <5(P + 9i - P"){P'X' p \J^\hA"(e z )){hA"(e z )\J v \PX p ) = (16) 

h A" 

= E'(^(P + P") £ (P'K\ ^|«A'(P'))(/ l A'(P')|/iA"(e z ))(/ l A"(e 2 )|/iA(P))(/ l A(P)| J„|PA P ). 

h A, A', A" 

Here \hA(a)) denotes the state with the spin projection onto vector a eq ual to A. 
Wave functions of these states are related via Wigner D-functions [llj : 

(hA"(e z )\hA(P)) =£fal(<p,0,0), (17) 

where ip, 8 are polar angles of the vector P and Sh is the spin of the state \h) . Using H17|) and properties of D-functions 
we have 

£>A'(i*)|M'V0><AA'^ (18) 

A" A" 

Using also amplitudes definition we obtain 

W»„{P'X p ;PX p ) = £ (-1) A+V 4 2 > A) x (19) 

A, A' 

xj:'(2n)±6(P + qi -P")f[ h \ q \) ^~M;f2A+l),v(2A' + i)(M,0), 

h 

where E\ and £2 are polarization vectors of the 1st (qi) and 2nd (92) photons of Eq. l|lU|) . defined according to (JSJ. 
After that the asymmetry becomes 



nD J 8s 

M 2 

Y^(^f5(P + k-P"- k") Y: f[ h \ql) }f{ql)X^{W, q Ul) 



A " = / ^r— dw2 / dn k „— x (20) 

irD J 8s J q fq% 

M 2 



A, A' 



where 

X^,(W,ll4) = * E (-^^^"^tlUAiPnutW-O) x (21) 

xi^ A) ef iV) fi A ,(P)(-7 S ^aK(i y ). 

The quantities X^, can be calculated explicitly (provided the proton form factors are known) , and the only unknown 
in 1)20(1 are electromagnetic transition amplitudes . 



III. THE MODEL FOR TRANSITION AMPLITUDES 



Obviously, it is practically impossible to take into account all allowed intermediate states. To proceed further, we 
need to restrict these states somehow. The authors of .9], for example, included only N and ttN states. Although 
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below the irirN threshold (£i a b ~ 0.3 GeV) such approach gives an exact result, one can expect that as the energy 
increases, this approximation becomes worse, since more intermediate states (e.g. rjN, iritirN) will be possible. 

In the present paper we use another way to model the intermediate states. We treat them as a number of resonances 
and neglect the non-resonant continuum contribution. At present we cannot estimate the non-resonant contribution 
well enough, but we can give qualitative arguments that it is small. 

At a glance one may conclude that the relative size of the non-resonant contribution will be approximately the 
same as in inelastic cross-sections or structure functions. But actually it is likely to be much smaller for the following 
reason. Contrary to strictly positive quantities, such as cross-sections, the asymmetry can have either sign. Thus the 
contributions from different non-resonant states will mostly cancel each other. This is similar to the fact that the 
average of many uncorrelated random quantities has much smaller dispersion that each of them. 

At P' = P the hadronic tensor W ^{P' \' p ] PX p ), which was introduced in previous section, turns into the hadronic 
tensor of inelastic ep-scattering. It is natural to assume that the qualitative properties of both tensors are similar, so 
we should first look what resonances contribute to the inelastic ep-scattering. 

There are three prominent resonant peaks in the inelastic ep cross-section: the so-called 1-st, 2-nd and 3-rd resonance 
regions. 

The first resonance peak is due to the A-resonance (P33(1232)), the second peak consists of Di3(1520) and Sn(1535). 
There are many resonances which contribute to the third resonance region, but there are serious arguments (see, e.g., 
Ref. that the dominant contribution comes from Fis(1680). Moreover, it is the only one, for which the transition 
amplitudes are known. Although the Roper resonance Pn(1440) does not contribute significantly to the inelastic 
ep-scattering 0, we also included it in our calculations. 

For the proton in the intermediate state one has 

]T (2tt) 3 6(P + qi - P") = J ^;S(P + qi - P") = S(W 2 - A/ 2 ), (22) 

h 

The resonance, however, has a mass Mr ^ M and some finite width Tr, so we "spread" (5-function with the relativistic 
Brcit-Wigncr formula: 

After that the expression '(2n) 3 S(P + q — P")/^(Si) /^(sl)) entering the formula for asymmetry, will take the 

h 

form 

The first part is the "elastic" (proton) contribution; in the second part the sum runs over all resonances taken into 
account. The quantities depend only on q 2 but not W. They are related to commonly used 01 ^3/2, ^i/2, an d 
S1/2 as 



/1 = kA 3/2 , i] R f-i = kA 1/2} and f = fjk ,,, ^1/2, (25) 

yWV - q HMr - M 2 - q 2 ) 2 



where n = \l m ^ AIr — ^— ^ , while for the proton they are 



7TU 



= 0, f { *\q 2 ) = 2MG E (q% f { _% 2 ) = -G M ( q 2 )^2f, (26) 

which is easy to derive comparing @ and 

For our calculation we use experimental data on Ah (i.e. j4 3 / 2 , ^1/2 : $1/2) given in |l3j| . Unfortunately, there are 
no data on Si/ 2 for Di3(1520) and i<i5(1680), so in our calculations we set it to zero. 

To evaluate the asymmetry we need to fit these data somehow. The fitting procedure is described in the next 
section. 
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IV. FITTING PROCEDURE 



Masses and widths of resonances were taken from PDG Q. The proton for m factors were modelled using the 
well-known dipole fit 

GM^)/, P = G E ^)= [l + ^ /Qir (27) 

where Q 2 = —q 2 , Qq = 0.71 GcV 2 and /x p ~ 2.79 is the proton magnetic moment. 

Now consider the fit of electroproduction amplitudes. According to quark model prediction the high Q 2 behaviour 
of the transition amplitudes should be like Ai/ 2 ~ 6*1/2 ~ Q~ 3 , A 3 / 2 ~ Q~ 5 US- 

For the proton we have 

^~^G u tf)~ + £ ~Q~\ (28) 



The denominator (dipole formula) is entirely due to quark structure, while the numerator is just a kincmatical factor 
1/ 



But the fact that A 1 / 2 tends to zero as Q 2 — > is the specific feature of the proton. For the other hadronic states 



Ai/ 2 (0) 7^ 0, so like [Tj| and for the same reasons we introduce factor yj (Mr — M) 2 + Q 2 and assume, instead of 

psli. 



y /(M R -M)* + * 

1/2 (i + Q 2 /Q 2 ) 2 • (29) 



So, to obtain correct asymptotic behaviour of Aw 2 at Q 2 — ► and Q — » 00, it is useful to fit the function 

which has finite values at both Q 2 = and Q 2 = 00. The same can be stated for Si/ 2 - On the other hand, A3/2 
has another asymptotic behaviour (~ Q~ 5 ), therefore we would have used (1 + Q 2 /Qq) 3 instead of (1 + Q 2 /Ql) 2 in 
expressions like (|29|) and (|30fl for A 3 / 2 . 

However the Q 2 values needed for our calculation are not too high. Since Q 2 = ^— — ^ sin 2 |, where 9 is the cm. 
scattering angle, one gets Q 2 lax ~ 3 GeV 2 at £?i a b = 2 GeV. Trying different parameterizations, we found that better 
agreement with the experimental data in the range Q 2 < 3 GeV 2 , especially for Delta, is achieved if we use the same 
formulae 



y/ (M R - M) 2 -f W 

{i + Q 2 /Ql? 



A h n 1 n 2/r>2\2 — ( 31 ) 



for all amplitudes. Thus we fit the functions 



To describe high Q 2 behaviour better, we treat Ah as a function of f = 1 — ,, n \ lr .- 2 instead of Q 2 . At low Q 2 it does 

not matter, since £ ~ Q 2 , but the advantage is that £ is finite at Q 2 — > 00, so we can use simple linear or polynomial 
fit for all £ values. 

In our calculations we restrict ourselves to the linear least-squares fit of the form 

A H (t)=a + b£. (33) 



Results of the fit are summarized in Table 1. The corresponding dependence of amplitudes Ah vs Q 2 for all considered 
resonances is shown on Figure together with experimental points. 
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State 


^4-3/2 -Al/2 S'i/2 


P 33 (1232) 
Dis(1520) 
Sn(1535) 
Fi 5 (1680) 
Pn(1440) 


-0.929+0.264 £ -0.485+0.130 £ 0.069+0.022 £ 
0.318-0.273 £ -0.029-0.474 £ no data 

0.123+0.416^ 0.212-0. 614c; 
0.185-0.052 £ -0.033-0.199 £ no data 

-0.351+0.787 £ 0.236-0.134^ 



TABLE I: Fit of the transition amplitudes. All values are in GeV 




FIG. 3: Fit of the transition amplitudes A 3 / 2 (solid lines and filled circles), ^1/2 (dashed lines and open circles) and S1/2 
(dash-dotted lines and crosses) for the resonancesfa) P 33 (1232), (b) 5n(1535), (c) Di 3 (1520), (d) Fi 5 (1680) and (e) Pn(1440). 
Experimental points are from a compilation of Hal . 
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V. NUMERICAL RESULTS AND CONCLUDING REMARKS 



Figure 0] displays the contribution of separate resonances to the target normal asymmetry A n vs the cm. scattering 
angle 9 at different electron lab. energy, Ey a ^. One sees that globally the A(1232) contribution is dominant. This 
is due to its large transition amplitudes, in comparison with other resonances, and the lowest mass among them. 
The contribution of the Roper resonance was obtained to be not negligible. Moreover at -Ei a b ~0.9 GeV it becomes 
comparable with the A(1232) contribution (the upper-right panel of Figure^}. This is very nontrivial fact, because 
the Roper contribution in inelastic eiV-scattering is very small. Nevertheless it can be studied in precise measurements 
of the A n at special kinematical conditions. 

One sees also that the contributions from the A and other resonances have mostly opposite sign and tend to cancel 
each other, especially at high beam energy. It is clearly seen from Figure where we plot the elastic (proton) and 
inelastic (resonance) parts of the asymmetry and the total asymmetry. The elastic contribution dominates at low 
energy (£i a b <0.3 GeV) and at energy higher than 1.3 GeV. It is quite obvious for the low energy, because the 
energy is insufficient for resonances to be produced. But at high energy it is nontrivial result, which has interesting 
consequences. As was discussed in Introduction, the asymmetry depends on the imaginary part of the amplitude. 
But since the real and imaginary parts are connected (via the dispersion relations), we may expect that the real part 
will also be defined mostly by proton contribution. This is important for the proper interpretation of the proton form 
factors measurements. 

In summary we have calculated the target normal spin asymmetry, A n , for the e~p — > e~p reaction at the electron 
beam energy up to few GeV in the laboratory frame. This quantity gives a direct information about the imaginary 
part of the reaction scattering amplitude and comes from the second and higher order pcrturbative terms. 

To calculate the imaginary part of the amplitude we used unitarity and saturated the intermediate hadron states by 
the proton (the so-called elastic contribution) and the resonances from the first, second and third resonance regions 
(the inelastic contribution). We neglect the non-resonant inelastic contribution, which we expect to be small (see 
section ITTl|) . Besides that, the calculated contributions of separate resonances are interesting alone. 

Our calculations demonstrate that at special kinematical conditions (the electron lab. energy near 0.9 GeV) 
the contribution of the Roper resonance Pn(1440) becomes comparable with the A(1232)-contribution and affects 
significantly the target asymmetry. It turn, this opens a possibility to study Pn(1440) electromagnetic transition 
amplitudes in precise measurement of the asymmetry. 
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FIG. 4: The contribution of resonances to the asymmetry at different electron lab. energies, (a) 0.57 GeV, (b) 0.855 GeV, (c) 
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1.4 GeV and (d) 2 GeV. Solid line — P 3 3(1232), long-dashed line 
line — Fi 5 (1680) and dotted line — Pn(1440). 
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FIG. 5: Target normal spin asymmetry for different electron lab. energies, (a) 0.57 GeV, (b) 0.855 GeV, (c) 1.4 GeV and (d) 
2 GeV. The dashed line is the elastic contribution, the dash-dotted line is the inelastic contribution, the solid line is total. 



